Constraining Bilinear i?-Parity Violation from Neutrino Masses 
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We confront the i?-parity violating MSSM model with the neutrino oscillation data. Investigating 
the 1-loop particle-sparticle diagrams with additional bilinear insertions on the external neutrino 
lines we construct the relevant contributions to the neutrino mass matrix. A comparison of the 
so-obtained matrices with the experimental ones assuming normal or inverted hierarchy and taking 
into account possible CP violating phases, allows to set constraints on the values of the bilinear 
coupling constants. A similar calculation is presented with the input from the Heidelberg-Moscow 
neutrinoless double beta decay experiment. We base our analysis on the renormalization group 
evolution of the MSSM parameters which are unified at the GUT scale. Using the obtained bounds 
we calculate the contributions to the Majorana neutrino transition magnetic moments. 
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SUPERSYMMETRIC MODEL WITH 
i?-PARITY VIOLATION 



The recent confirmation of neutrino oscillations [l| 
gives a clear signal of existence of physics beyond 
the standard model of particles and interactions (SM). 
Among many exotic proposals the introduction of su- 
persymmetry (SUSY) proved to be both elegant and 
effective in solving some of the drawbacks of the SM. 
The minimal supersymmetric standard model (MSSM) 
(a comprehensive review can be found in Q) populates 
the so-called desert between the electroweak and the 
Planck scale with new heavy SUSY particles, thus remov- 
ing the scale problem. What is more, using the MSSM 
renormalization group equations for gauge couplings in- 
dicates that there is a unification of g\ , gi and (73 around 
tugut ~ 1.2 x 10 16 GeV which means that MSSM in 
a somehow natural way includes Grand Unified Theories 
(GUTs). This model is also characterized by a heavier 
Higgs boson, comparing with the Higgs boson predicted 
by SM, which is in better agreement with the known ex- 
perimental data. New interactions present in MSSM lead 
to many exotic processes which opens a completely new 
field of research. 

Building the minimal supersymmetric version of the 
Standard Model one usually assumes the conservation of 
the i?-parity, defined as R — (— i) 3 - B + i + 2S j where B is 
the baryon number, L the lepton number, and S the spin 
of the particle. The definition implies that all ordinary 
SM particles have R = +1 and all their superpartners 
have R = — 1. In theories preserving i?-parity the prod- 
uct of R of all the interacting particles in a vertex of 
a Feynman diagram must be equal to 1. This implies 
that the lepton and baryon numbers are conserved, and 



that SUSY particles are not allowed to decay to non- 
SUSY ones. It follows that the lightest SUSY particle 
(usually the lightest neutralino Xi) must remain stable, 
giving a good natural candidate for cold dark matter. All 
this makes the i?-parity conserving models very popular. 

In practice, however, the i?-parity conservation is 
achieved by neglecting certain theoretically allowed terms 
in the superpotential. Casting such hand-waving ap- 
proach away, one should retain these terms, finishing 
with an i?-parity violating (RpV) model, with richer phe- 
nomenology and many even more exotic interactions 0- 
6]. The RpV models provide mechanisms of generat- 
ing Majorana neutrino masses and magnetic moments, 
describe neutrino decays, SUSY particles decays, exotic 
nuclear processes like the neutrinoless double beta decay, 
and many more. Being theoretically allowed, RpV SUSY 
theories are interesting tools for studying the physics be- 
yond the Standard Model. The many never-observed pro- 
cesses allow also to find severe constraints on the non- 
standard parameters of these models, giving an insight 
into physics beyond the SM. 

The violation of the i?-parity may be introduced in 
a few different ways. In the first one i?-parity violation 
is introduced as a spontaneous process triggered by a 
non-zero vacuum expectation value of some scalar field 
Another possibilities include the introduction of ad- 
ditional bi- [I| or trilinear @ RpV terms in the super- 
potential, or both. In the following we incorporate the 
explicit RpV breaking scenario. 

The i?-parity conserving part of the superpotential of 
MSSM is usually written as 
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while its RpV part reads 
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and M^o is the standard MSSM neutralino mass matrix: 



The Y's are 3x3 Yukawa matrices. L and Q are the 
SU(2) left-handed doublets while E, U and D denote the 
right-handed lepton, up-quark and down-quark SU{2) 
singlets, respectively. H d and H u mean two Higgs dou- 
blets. We have introduced color indices x,y,z — 1,2,3, 
generation indices k = 1, 2, 3 = e, fj,, r and the SU(2) 
spinor indices a, b = 1,2. 

As far as the (in principle unknown) RpV coupling 
constants are concerned, the most popular approach is 
to neglect the bilinear terms and to discuss the effects 
connected with the trilinear terms only. In such a case, 
if one is not intereseted in exotic baryon number violat- 
ing processes, one has to additionally set A" = 0, which 
ensures the stability of the proton. In this paper we con- 
centrate on the bilinear terms only and set all trilinear 
RpV couplings to zero. 

For completeness we write down the scalar mass term 
present in our model, 

C nass = m ^/ 1 t/ ld + m ^^ u + g t m | g + / t m 2; 

+ uva^jv) + dva^d) + em E e f , (3) 
the soft gauginos mass term (a = 1, 8 for gluinos) 

C gaug. = (m^b + M 2 W i *W i + M 3 g~Jg a + h.c. 

as well as the supergravity mechanism of supersymmetry 
breaking, by introducing the Lagrangian 
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+ (Au) ijq rh b u u jx + Bfih a d h b u + B lKl l?h b u ,(5) 



where lowercase letters stand for scalar components of the 
respective chiral superfields, and 3x3 matrices A as well 
as Bfi and Bi are the soft breaking coupling constants. 
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The matrix ^ has the seesaw-like structure and con- 
tains the sneutrino vacuum expectation values (vevs) Wj. 
These are in general free parameters which contribute to 
the gauge boson masses via the relation 



E 

i— e,/i,r 



2 2 
OJj = V = 



(246 GeV) 2 



(9) 

where vi and w 2 are the usual down-type and up-type 
Higgs boson vevs, respectively. By introducing the angle 
/3 defined by tan /3 = v 2 jv\ we obtain four free param- 
eters of the theory: tan/3 and u>i. Fortunately it turns 
out that in order to obtain proper electroweak symmetry 
breaking the sneutrino vevs cannot be arbitrary. We give 
the details in the next section. 



III. HANDLING THE FREE PARAMETERS 

The RpV MSSM model introduces several new free 
parameters when compared with the usual MSSM. For- 
tunately their number can be constrained by imposing 
GUT unification and renormalization group evolution. 
In this paper we restrict ourselves to the bilinear RpV 
couplings only, setting all trilinear couplings (A, A', A") 
to zero. This assumption simplifies some of the RGE 
equations, which we list below. Such approach leads at 
the end to the following set of free parameters: too, to^, 
A , tan/3, sgn(^i), and nf UT (i = 1,2,3). 



A. Masses and soft breaking couplings 



II. NEUTRINO— NEUTRALINO MIXING 

The inclusion of the bilinear RpV terms imply mix- 
ing between neutrinos and neutralinos. In the ba- 
sis (v e ,Ufj,, v T ,B, W 3 ,H^,H°) the full 7x7 neutrino- 
neutralino mixing matrix may be written Q in the fol- 
lowing form: 

where 

(-ig'w e \gw e -K e \ 

-\g'up. \gu„, -Kn (7) 

-\g'u T \guj T -n T J 



The masses of all the supersymmetric scalars are uni- 
fied at ttigut to a common value mo, and of all the super- 
symmetric fermions to to^- The values of the trilinear 
soft SUSY breaking couplings are set according to the 
following relations |8| 

Ae,d,u — A Ye,d,u, (10) 
B = B 1>2>3 = A -l. (11) 

The RGE equations for the A couplings can be found 
elsewhere pHla]. The B couplings are evolved down to 
the low energy regime according to the renormalization 
group equations 

16n2 ll = 6 Tr(A[/Y ^ + 6 T < AdY d) 

+ 2 Tr(A B Y^) + &g 2 2 M 2 + 2g 2 1 M 1 , (12) 
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FIG. 1: An example of RG running of the bilinear Ki couplings. The unification scenario was: mo = 200 GeV, m 1( / 2 = 500 GeV, 



A a = 200, sgn(/x) = 1. All 



,GUT 



were equal to 1 MeV. 



16tt 



6 Tr(A[/Y[,) + 6<? 2 M 2 + 2&Mi, (13) 



dt 

= 6 Tr(A[/Y^) 

+ 2 Tr(A E V\g) + §9\M 2 + 1g\M x . (14) 



where g\ = 5/3 g' 2 /(47r 2 ) and g 2 = g 2 /(4:ir 2 ), 5/3 being 
the GUT normalization factor. 



B. Bilinear Ki couplings 

The three K,f UT couplings at GUT scale remain free in 
our model. After setting them the couplings are evolved 
down to the mz scale according to the renormalization 
group equations which in our case take the following sim- 
ple form: 
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k 4 (3 Tr(YuYl 

3 



352 



al) 



(15) 



An example of the running of m is presented on Fig. [T] 
One sees that for higher tan /3 the couplings vary rather 
weakly (notice the logarithmic scale on the energy axis) 
for the whole energy range between the electroweak scale 
mz and itigut- For small tan/3 < 10 the difference 
between the itigut and mz values are of the order of 
< 35%. The value 1 MeV at the GUT scale was chosen 
arbitrarily; we will show later that this is the typical or- 
der of magnitude for which agreement with experimental 
data on neutrino masses and mixing may be obtained. 



C. Vacuum expectation values 

At the beginning of the numerical procedure we set the 
down and up Higgs vevs to 

vi = v cos f3, v 2 = v sin /?, (16) 

while the initial guess for the sneutrinos vevs is 

Ui =0. (17) 

The actual values of 0Ji are calculated from the condition 
that at the electroweak symmetry breaking scale the lin- 
ear potential is minimized. By taking partial derivatives 
of the potential one obtains the so-called tadpole equa- 
tions [8], which are zero at the minimum. 

In our procedure we solve three equations, which can 
be written as (i = 1, 2, 3) 



3=1 



where 



Qji = KjKi + (m|)ji + SjiD, 



Sji being the Kronecker delta, and 

D= 1 -(g 2 +g' 2 )(v 2 -2vl). 



(18) 



(19) 



(20) 



Notice that they are linear in 0^2,3 and therefore this 
set has only one solution. After finding it, we use the 
trigonometric parameterization Q, which preserves the 
definition of tan/3, 



vi 



v sin a± sin a 2 sin 0:3 cos /3, 



(21) 
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FIG. 2: (a) The basic 1-loop diagram giving rise to the Ma- 
jorana neutrino mass in the i?-parity violating MSSM. (b) 
1-loop diagram with RpV neutrino-neutralino couplings in- 
cluded on the external lines. 





FIG. 3: Diagrams with bilinear neutrino-neutralino interac- 
tions leading to the Majorana neutrino mass. 



i)2 = vsinai sina^ sin 013 sin/3, (22) 

UJ\ = wcosai sina2 sina3, (23) 

0J2 = wcosa2sina3, (24) 

0J3 = v cos 0:3, (25) 

to calculate new values of vi and ■ We return back 
to the tadpoles with these new values and continue in 
this way until self-consistency of the results is reached. 
It turns out that due to the expected smallness of the uii 
vevs, the initial guess Eq. (fTT|) is quite a good approx- 
imation. It usually suffices to repeat the whole proce- 
dure three times to obtain self-consistency at the level of 
0(1O -4 ), which is more than enough for our purposes. 
The so-obtained set of vevs is used during the determi- 
nation of the mass spectrum of the model. 



IV. FEYNMAN DIAGRAMS WITH RPV 
COUPLINGS ON THE EXTERNAL NEUTRINO 
LINES 



TABLE I: Nine diagrams with neutrino-neutralino mixing on 
the external lines leading to Majorana neutrino mass. / and 
III are the neutralinos which mix with the neutrinos on the 
external lines of the diagram depicted on Fig. [3] // defines 
the content of the loop (uu stands for the up-type (s)quarks, 
dd stands for the down-type (s)quarks,gg for all (s)quarks, 
and ll for (s)leptons). 
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It is well known that, once allowing for i?-parity viola- 
tion, a particle-sparticle 1-loop diagrams give important 
corrections to the usual tree level neutrino mass term. 
These processes have been extensively discussed in the 
literature [l6l - |22l |. mainly in the context of constraining 
the tree-level alignment parameters A or the trilinear 
RpV couplings A and A' 0,111. 

In general, the explicit RpV effects may be taken into 
account in three different ways. One may include the bi- 
linear RpV couplings or the trilinear couplings, or both. 
Of course the most complete one is the third possibility, 
which is at the same time the most complicated. There- 
fore it is customary to limit the discussion to either tri- 
or bilinear terms only. In this paper we are interested in 
bilinear couplings and set all trilinear couplings to zero. 

In order to discuss the possible magnitude of the bilin- 
ear RpV couplings Hi we extend the simplest diagrams by 
including the neutrino-neutralino mixing on the external 
lines. 

The topology of the basic type of 1-loop diagrams 
we will consider is presented on Fig. [UJa). These dia- 
grams lead to Majorana neutrino mass term, where the 
effective interaction vertex is expanded into the RpV 
particle-sparticle loop. These diagrams and their more 
complicated versions with the Higgs bosons and sneutri- 



nos inside the loop, were classified in e.g. Ref. [19| and 
discussed in details elsewhere (see [16l - f23j among oth- 
ers). In the present paper we add the possible neutrino- 
neutralino mixing on the external lines (Fig.[2jb)), which 
leads to another contributions to the neutrino mass. Ob- 
viously this additional contribution must be in agreement 
with the present experimental data. We discuss two main 
cases, in which either lepton and slepton or quark and 
squark are in the loop (in the case of higgsino H u only 
the up- type quarks count). At the same time the neu- 
trino may mix either with the gauginos: bino B or wino 
W 3 , or with the neutral up- type higgsino H u . All the 
nine cases together with the relevant bi- and trilinear 
coupling constants have been gathered in Tab. HI 

The contributions from individual diagrams have been 
calculated using the same technique as in Refs. [20I [2lj|. 
In Ref. [2l| we have discussed the possible influence of 
including the quark mixing in the calculations. Here we 
neglect this effect. 

The neutrino mass matrix resulting from the bilinear 
processes only can be written as the following sum: 

9 

M ab = J2Kb> (26) 

i=l 
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where the separate contributions read 
1 C1C2C3C4 



M 



ab 



167T 2 mjmiu 



(27) 



The masses of the neutralinos mj and mjjj, and the cou- 
pling constants have to be taken from Tab. |TJ The func- 
tions F represent the contributions from the particlc- 
sparticle loops. They read: 



F, 



E 

F<ld = E 



3 8m 2 26 ' m u i f (x% ,x\ J ) 

„ sin 20-* y jj 

m d 'J{x 2 f,X{ J ) 

m q if(x% ,x l -l) 



2 

sin20 J ' 
sin 2<^ 



(28) 
(29) 
(30) 
(31) 



where and ^ are the j-th squark and slepton mass 
eigenstates' mixing angles, respectively. For simplic- 
ity we have defined dimensionless quantities ij^ = 
(m q a/m^b 2 ) 2 , which are the et-th quark mass over the 
6-th squark first or second mass eigenstate ratios. An 



analogous expression involving the lepton and slepton 
masses has been named yf b . The function coming from 
integrating over loop momentum is f(x, y) = [log(y)/(y— 
1) — log(a;)/(x — 1)]. The j-sums run over all squarks in 
F U ui Fddi an< ^ Fqq, and over all sleptons in F,r. The i- 
sums count all quarks in F qq - , up- type quarks only in F u a , 
down-type quarks in Fjs, and all leptons in F,r. The fac- 
tor 3 in F U u, J^jj, and F qq accounts for summation over 
quarks' colors. It is absent in the case of leptons. 

We do not discuss the M 1 contributions separately. 
The reason is that for different cases the couplings C2 
and C3 enter with opposite signs causing cancellations 
between such terms. Since none of the M % can show 
up without the others, only the full sum Eq. ([26| gives 
a meaningful picture. 



PHENOMENOLOGICAL MAJORANA 
NEUTRINO MASS MATRIX 



The neutrino mass matrix can be constructed from 
the Pontecorvo-Maki-Nakagawa-Sakata mixing matrix 
Upmns under certain assumptions. The matrix Upmns 
is usually parameterized by three angles and three (in the 
case of Majorana neutrinos) phases as follows: 
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(32) 



where s ?1 = sin( 



i 1 Ci 



cosdij. Three mixing angles 9ij (i < j) vary between and tt/2. The S is the CP violating 



Dirac phase and <p2, (j>3 are CP violating Majorana phases. Their values vary between and 27r. The explicit expression 
for the phenomenological mass matrix M v \ in terms of m t , dij, S, 4>2, 4>3 is given by 21]: 
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M TT = (s?, 3 sl 2 - 2 C 2 3 C12 S 23 S12 S U e 1 + <?23<?xi S \z { 
+ ( c 12 s 23 + 2 C 23 C12 S23 Sl2 Sl3 e^" 5 + sf 2 sf 3 < 



-2i8\ 



mi 

m 2 e- 12 ^ + 4 3 cj 3 m 3 e~ i2 ^ 



(33) 



In order to calculate numerical values of elements of 
this matrix one needs some additional relations among 
the mass eigenstates 7711,2.3- Experiments in which neu- 
trino oscillations are investigated allow to measure the 
absolute values of differences of neutrino masses squared 
and the values of the mixin g an gles. The best-fit values 
of these parameters read [H |24| 



\m x -m 2 \ = 
1 2 2 

sm 2 {6 12 ) = 

sin 2 (0 23 ) = 



713 



7.1 x 10~ 5 eV 2 , 

2.1 x 1(T 3 eV 2 , 

0.2857, 

0.5, 

0. 



(34) 



The present experimental outcomes are in agreement 
with two scenarios: 

• the normal hierarchy (NH) of masses imply the re- 
lation 777 1 < 777-2 < 7773, 

• the inverted hierarchy (IH) of masses imply the re- 
lation 7773 < 777l < 7772- 

Notice that in order to keep the same notation for the 
differences of masses squared and the mixing angles, the 
neutrino mass eigenstates are labeled differently in the 
NH and IH cases. 

At this point we are left with four undetermined pa- 
rameters, which are the phases and the mass of the light- 
est neutrino. To obtain most stringent limits on the new 
physics parameters the later is taken to be zero. As far as 
the phases are concerned we consider two separate cases. 
First we take all possible combinations of phases and for 
each entry of the matrix we pick up its highest possible 
value. In this way we obtain unphysical matrices, which 
give however some idea about the upper limits on the 
non-standard parameters. The maximal matrices for the 
NH and IH scenarios read as follows [2l|: 



K 'I max — 



0.00452 
0.00989 
,0.00989 



|X|(IH) < 



0.0452 
0.0312 
,0.0312 



0.00989 
0.02540 
0.02540 



0.0312 
0.0240 
0.0239 



0.00989 \ 

0.02540 eV, (35) 
0.02540/ 



0.0312X 

0.0239 

0.0240/ 



eV. 



(36) 



The more conservative approach assumes that the CP 
symmetry is preserved which can be achieved by neglect- 
ing the phases present in the Upmns matrix. In such 
a case the NH and IH matrices take the following forms: 



\M 



(NH) _ 



'0.00240 0.00269 0.00269 \ 
0.00269 0.02553 0.01951 eV, (37) 
,0.00269 0.01951 0.02553/ 



\M 



(IH) 



'0.045267 
0.000249 
,0.000249 



0.000249 
0.022801 
0.022801 



0.000249 \ 
0.022801 eV. 
0.022801/ 



(38) 

Yet another possibility is to construct M. using con- 
straints from non-observability of the neutrinoless double 
beta decay (0^2/3). The study of the 0^2/3 decay [H[ is 
one of the most sensitive ways known to probe the abso- 
lute values of neutrino masses and the type of the spec- 
trum. The most stringent lower bound on the half-life 
of 0^2/3 decay were obtained in the Heidelberg-Moscow 



7G 



Ge experiment (H (T®J~ exp > 1.9 x 10 25 yr). By as- 
suming the nuclear matrix element of Ref. [27| we end up 
with |t77^| = U 2 X mi + U 2 2 ni2 + U 2 3 7773 < 0.55 eV, where 
U is the neutrino mixing matrix Eq. (|32p . The element 
\mpp\ coincides with the ee element of the neutrino mass 
matrix in the flavor basis and fixing it allows to construct 
the full maximal matrix, which reads: 



\M 



(HM 




1.29 
1.35 
1.04 



eV. 



(39) 



In the next section we present the results for each of these 
five cases. 



VI. CONSTRAINING k COUPLINGS FROM 
THE NEUTRINO MASS MATRIX 

Our aim is to find constraints on the Hi coupling con- 
stants coming from the neutrino mass matrices. As an 
example of the unification conditions we take the follow- 
ing input: 

A a = 200, m = 200 GcV, 777i /2 = 500 GeV, (40) 
and additionally 



tan/3 = 10, sgn(/j) = 1. 



(41) 



We do not expect great differences in the results if the 
GUT conditions were changed. The only exception may 
be the tan /3 parameter (defined at 777 z scale) which dom- 
inates the running of the k's. By looking on Fig. [1] only 
very low values of this parameter will influence the results 
significantly. 
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TABLE II: Some results for the SUSY scenario tan /9 = 10, A = 200, m = 200 GeV, m 1/2 = 500 GeV. 



GUT GUT GUT 
Kx K 2 K 3 


Resulting mass matrix 


Compare with 


Remarks 


[MeV] 


[eV] 







9.50 14.80 14.; 



9.46 13.02 13.02 



0.85 2.03 2.03 



2.72 1.98 1.' 



0.62 2.03 2.03 



0.62 0.70 0.70 



2.72 1.92 1.93 



0.27 0.19 0.19 



0.553927 
0.865038 
0.859783 

* 0.548762 
0.757853 
0.752398 

* 0.004520 
0.010734 
0.010686 

* 0.045316 
0.032945 
0.032951 

* 0.002402 
0.007821 
0.007787 

* 0.002402 
0.002688 
0.002684 

* 0.045316 
0.032065 
0.032084 

* 0.000453 
0.000320 
0.000320 



864372 
349844 
341638 
757606 
046272 
038740 
010728 
025474 
025361 
032954 
023958 
023963 
007824 
025474 
025361 
002691 
003011 
003007 
032053 
022680 
022694 
000320 
000226 
000226 



0.861090 

1.344714 

1.336555 

0.753587 < 

1.040721 

1.033237 

0.010698 

0.025404 

0.025292 

0.032976 { 

0.023974 

0.023978 

0.007802 ( 

0.025404 

0.025292 

0.002687 < 

0.003007 

0.003003 

0.032118 ( 

0.022726 

0.022740 

0.000321' 

0.000227 

0.000227 



max 

p,T elements to big 



M 



(HM) 
max 



M 



(NH) 



M (m) 



M 



(NH) 



efj, and er elements to big 
efj, and er elements to big 
efi, er and fir elements to big 

X(NH) 

A4^ IH ^ e/x and er elements to big 






I I I I I I I 




• •- '^%--'"V"tv--- 


: ' )* 


>""y* <**•-'■ ''■>' '■ .'h'*'.^i *■ 






■t" ! '..' 


;>%uf- '* ^ *.**'*■•_-, '•. J 


v"/ ■■■ 




3 4 6 8 1 
Kj IMeV] 



1 *&%■■*>■■ 

4 fj^^^^e 




i , 2 l- ; ',r-. *'/'c 



0.2 0.4 



0.4 0.8 1.2 1.6 2 

k 2 [MeV] 





FIG. 4: Allowed parameter space in the maximal HM case. 



FIG. 5: Allowed parameter space in the maximal NH case. 



We proceed in two steps. Firstly we find such values of 
nf UT which will reproduce the diagonal elements of the 
mass matrices. This can be achieved with good accuracy, 
but it turns out that some of the elements (off-diagonal, 
see remarks in Tab. |TTJ) of the resulting matrix exceed the 



allowed values. It means that the n's will not take their 
maximal values simultaneously. 

Secondly we go down with the K,f UT to lower the off- 
diagonal elements to the acceptable level. This, however, 
can be done in many different ways. Some explicit ex- 
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05 1 15 3 04 08 1.2 • fi 2 



FIG. 6: Allowed parameter space in the maximal IH case. 




04 08 1 2 1.6 2 



FIG. 7: Allowed parameter space in the NH case with con- 
served CP symmetry. 

amples are listed in Tab. [TT| but to find the full allowed 
parameter space we have prepared scatter plots which 
are presented on Figs. 0HS1 Each of the plots consists of 
roughly 2000 points chosen randomly from the intervals 
between zero and 1.1 times the assessed upper limit for 
given K,f UT . 

The boundaries of the allowed parameter space for 
Ki in the case of unphysical neutrino mass matrices 
|A / f|nmx I ' NH ' IH ' > are nearly box-shaped, except for a small 
region of excluded values in the upper right-hand corner 



FIG. 8: Allowed parameter space in the IH case with con- 
served CP symmetry. 

of some of the plots. This behavior is expected, as the 
mass matrices to be reproduced contain on each entry 
the maximal allowed value for it. 

Much more interesting shapes are obtained for the CP 
conserved cases (Figs.[7]and|8]). The projections onto the 
(«i, K2) and (ki, K3) planes are nearly identical and con- 
tain non-linear boundary parts. The most constrained is 
the inverted hierarchy case with conserved CP (Fig. [8]) 
due to two orders of magnitude differences between the 
diagonal fifj, and rr and off-diagonal e\i and er elements. 



VII. THE TRANSITION MAGNETIC MOMENT 

The RpV loop diagrams provide not only an elegant 
mechanism of generating Majorana neutrino mass terms, 
but also, after a minor modification, may be the source 
of the transition magnetic moment \j, a \,. This quantity 
represents roughly the strength of the electromagnetic 
interaction of the neutrino. Since the latter is electri- 




FIG. 9: Feynman diagram with neutrino-neutralino mixing on 
the external lines, leading to the Majorana neutrino transition 
magnetic moment. 
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TABLE III: Contribution to the Majorana neutrino transi- 
tion magnetic moments coming from the bilinear neutrino- 
neutralino mixing, for the GUT scenario: Aq — 200, mo = 
200 GeV, mi/a = 500 GeV, tan/3 = 10. 



/'{: 7 



trilinear only 



-20 
20 



IH-CP 7.0 x 10 -22 7.0 x 10 -22 6.0 x 10 

IH-max 8.8 x 10~ 20 8.5 x 10~ 20 6.5 x 10 

NH-CP 7.6 x 10~ 21 7.6 x 10~ 21 5.5 x 10 

NH-max 2.8 x 10~ 20 2.8 x 10~ 20 7.0 x 10 

HM-max 2.4 x 10~ 18 2.3 x 10" 18 2.9 x IP" 18 < 10 



-20 
-20 



< 10 

< 10 



< io 

< 10 



-19 

-17 



-18 
-17 



-15 



cally neutral, the interaction must take place between an 
external photon and a charged particle from inside the 
virtual RpV loop. In practice, only the photon-fermion 
interactions are taken into account, since the photon- 
boson (squark or slepton) interaction would be strongly 
suppressed by the big mass of the SUSY particle. The 
relevant Feynman diagram is presented on Fig. |9] 

The contribution to the Majorana neutrino magnetic 
moment from the discussed diagrams is given by (in Bohr 
magnetons hb) 



Hab 



(1 



id 



Jab) 



m„i 



4tt 2 

Qq* 



n C-2.Cz 



mimin 



w 



Ms- 



(42) 



Here we have denoted the electric charge of a particle (in 
units of e) by Q. The dimensionless loop functions w 
take the forms 

( g ) sm26» J y ij (£) sin20 J y y 
w ij = — ^ — 3{x 2 ,x{), w\/ = — - — g{y^,y{), 

(43) 



where 6, <fi, X1.2, and y\^, are the same as in Eqs. 
(l3Tj) . and g(x, y) = (x log(x) - x + 1)(1 - x)~ 2 - (x -4 y). 
The sum over i and j in Eq. (|42[) accounts for all the 
possible quark-squark and lepton-slepton configurations 
for given neutralinos. The factor 3 in front of counts 
the three quark colors. 

The results for the already discussed GUT parameters 
are presented in Tab. IIIIl The last column contains for 
comparison upper bounds for the magnetic moment in 
the case when only trilinear interactions are taken into 
account. One sees that they are at least one order of 
magnitude stronger than the discussed bilinear contribu- 
tions. In Tabs. ED and [V] we show the results of similar 
calculations for two other set of parameters. In Tab. IIVI 
the unification parameters are 'low', while in Tab.fVlthcir 
values are 'higher'. The last column is given as previously 
for comparison. Also here the conclusion is clear, that the 
discussed contribution to the main process is at best of 
the same order of magnitude, in most cases being at least 



TABLE IV: Like in Tab. [TTT] but with Ao = 100, m = 
150 GeV, m 1/2 = 150 GeV, tan/3 = 19. 



Me 



/'<: 



/i MT trilinear only 



IH-CP 3.0 x 10" 21 2.9 x 10~ 21 2.5 x 10~ 19 < 10~ 18 

IH-max 3.7 x 10" 19 3.6 x 10" 19 2.7 x 10~ 19 < 10" 18 

NH-CP 3.2 x 10" 20 3.1 x 10~ 20 2.2 x 10~ 19 < 10" 18 

NH-max 1.2 x 10" 19 1.1 x 10~ 19 2.9 x 10~ 19 < 10" 18 

HM-max 1.0 x 10~ 17 9.8 x 10~ 18 1.2 x 10~ 17 < 10~ 16 



TABLE V: Like in Tab. |HJ] but with A = 500, m 
1000 GeV, mi/a = 1000 GeV, tan ,5 = 19. 



He 



[i^ T trilinear only 



IH-CP 3.7 x 10~ 22 3.7 x 10 -22 3.3 x 10 _2 ° < 10 _2 ° 

IH-max 4.6 x 10" 20 4.6 x 10~ 20 3.5 x 10~ 20 < 10~ 20 

NH-CP 4.0 x 10" 21 4.0 x 10" 21 2.9 x 10" 21 < 10" 20 

NH-max 1.4 x 10~ 20 1.5 x 10~ 20 3.7 x 10" 20 < 10~ 20 

HM-max 1.3 x 10 -18 1.2 x 10 -18 1.5 x 10" 18 < 10 -18 



an order of magnitude weaker. The reason for such a situ- 
ation is due to the high masses of the neutralinos, which 
enter the formula Eq. (l42t in the denominator. It was 
possible that they will be compensated by the unknown 
coupling constants proportional to the sneutrino vacuum 
expectation values u, especially that some of them are 
multiplied by negative numbers, cf. Tab. [TTJ Our explicit 
calculation showed that it is not the case. Also the ob- 
served differences between the values of fi a b, reaching not 
more than one order of magnitude, are mainly due to the 
changed value of the parameter tan /?, and only partially 
due to different values of the remaining parameters. 



VIII. SUMMARY 

The -R-parity violating MSSM has many free parame- 
ters which lower its predictive power. On the other hand 
this fact makes the model very flexible. In this paper 
we have presented a method of constraining the bilinear 
RpV couplings n. 

We have calculated the contributions to the neutrino 
mass matrix coming from the neutrino-neutralino mixing 
in processes in which the effective vertex is expanded into 
a virtual quark-squark or lepton-slepton loop. These 
contributions have been compared with the phenomeno- 
logical mass matrices derived using the best-fit experi- 
mental values of the neutrino mixing angles and differ- 
ences of masses squared. We discuss four cases in which 
normal and inverted hierarchy is explored both with con- 
served CP symmetry and with maximal values of each 
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matrix element. We also present the fifth case in which 
the neutrino mass matrix is calculated from the data 
published by the Heidelberg-Moscow neutrinoless dou- 
ble beta decay experiment. 

In general we have found that setting the k couplings at 
the unification scale to values of the order of < 0(1 MeV) 
renders the mass contributions correctly below the exper- 
imental upper bound. Another observation is that the 
bilinear RpV mechanism alone is not sufficient to repro- 
duce the whole mass matrix. This is, however, acceptable 
because in the general RpV loop mechanism one has to 
sum up the contributions from the tree-level (l6j . 

Mi + M 2 tan 2 6 W 

X 4(/*m^(Mi + M 2 tan 2 9 W ) sin 2/3 - M X M 2 ^) ' 

(44) 

where Aj = (iuji~VdKi are the so-called alignment param- 
eters, as well as contributions coming from the 1-loop 
diagrams (see Fig. Eta)), which are proportional to the 
totally unconstrained trilinear couplings A and A'. These 
parameters may be easily fine-tuned to reproduce the full 
mass matrix and we shift this discussion to an upcoming 
paper. 

The knowledge of the bounds on the K coupling con- 



stants allows one to discuss many exotic processes, like 
the neutrino decay and the interaction of neutrino with 
a photon, to mention only a few. The former may oc- 
cur as a two-step process, first through bilinear mixing 
with neutralinos, and then the decay of the actual neu- 
tralino. The later has been presented in the previous 
section showing by explicit calculation that this contri- 
bution does not exceed the main 1-loop mechanism. 

In our calculations we have fixed the GUT unification 
parameters. Due to technical difficulties in performing 
a full skan over the allowed parameter space we have 
picked only three representatives for which the calcula- 
tions were performed. We expect that the results will 
not change qualitatively with the changes of the input 
parameters, which is of course an assumption that may 
be worth checking. 
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